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ABSTRACT 


In this report we investigated the variation of dose rate with solid 
angle fraction in structures of simple geometries and compared the ex¬ 
perimentally determined infinite field dose rate with that calculated 
using the formula 

D* = B 0. 5 S E |S (u )-S (u ,)| + 0.088(1-S ) 
w w 1 a' u 7 a' u 7 J v w 7 

The geometry factor for skyshine radiation S (to) was taken from Figure 
B37 of "Structure Shielding Against Fallout Radiation from Nuclear 
Weapons" by L. V. Spencer and was applied to wall-scattered radiation 
as described in "An Engineering Method for Calculating Protection Af¬ 
forded by Structures Against Fallout Radiation" by C. Eisenhauer. The 
factors B w (X e ), E, and S w (X e ) were taken from the Engineering Manual 
"Design and Review of Structures for Protection from Fallout Gamma 
Radiation." 

Except for the thinnest walls, the values of experimental and calcu¬ 
lated dose rates were within 20%. This was greater than the estimated 
experimental error, and the difference followed a pattern that suggests 
a definite discrepancy between experiment and calculation. The same 
pattern of discrepancies was found in a direct comparison of experi¬ 
mental and theoretical values of S (w). 
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S (w ) - S (a) ,) 
a v u 7 a' u 7 j 
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CHAPTER 1 


INTRODUCTION 

The prediction of dose rates to be expected in a radiation shelter due to exter¬ 
nal sources presents a problem far more complicated than that of determining the 
attenuation of a narrow, collimated beam of radiation traversing a thickness of 
material and being measured by a collimated detector. The dose rate in a shelter, 
besides the narrow beam contribution of direct radiation, will include secondaries 
scattered through all angles and often of comparable intensity. 

During recent years this aspect of radiation shielding has been intensively 
investigated with particular emphasis on applications to shelters. Of particular 
importance has been a series of machine computations of gamma ray scattering and 
attenuation functions prepared and compiled by Dr. L. V. Spencer 1 using the moment 
method and Monte Carlo techniques based on fundamental cross-section data. Com¬ 
putations of the angular distribution of scattered radiation at various distances from 
point and plane sources yielded twenty-five graphs from which flexible, self- 
consistent methods of analysis were derived by Mr. Charles Eisenhauer and Mr. 

L. N. FitzSimons working in conjunction with Dr. Spencer. These curves, com¬ 
puted for Co-60, Cs-137, and 1.12-hr fallout radiation, have been published in the 

2 

so-called Engineering Manual , which is designed to assist architects and engineers 
in the design of structures for fallout protection. 

To date, these methods have been experimentally evaluated only by whole building 
exposures — structures so complicated as to yield little basic information, though 
providing significant confirmation of the general applicability of the method. The 
basic premises of this important field, therefore, have not been tested experimen¬ 
tally. This experiment was planned to test the premises by using the simplest possi¬ 
ble structure, thus reducing the measured dose rate to a function solely of the sub¬ 
tended solid angle and the thickness of the scattering walls. 

A cylindrical shell at the center of a circular field was used to provide this sim¬ 
plicity, with the additional advantage of circular symmetry, making possible the 
replacement of the plane source by a succession of point sources along one radius. 
Wall thickness was changed by adding successive half-cylinders, and solid angle 
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was changed by placing the detector at different depths in a hole beneath the cylinder. 
Below-surface positions were necessary to eliminate direct radiation. The infor¬ 
mation derived from this experiment about the variation of dose rate with solid angle 
for a given cylinder thickness could then be compared with calculations based on the 
Engineering Manual and Spencer’s Monograph. In the second part of the experiment, 
the cylinder was replaced by a wall of 1-in. thick steel, 6 ft by 12 ft, to test the 
possible influence of structure configuration on the dose rate. 

The experiment is described in detail in Chapter 2. The basic data are pre¬ 
sented in Chapter 3 and are analyzed in Chapter 4. Chapter 5 contains the conclu¬ 
sions derived from the experiment. An experimental curve for S (cu) is derived in 

cl 

Chapter 6 and is compared with the theory. 







CHAPTER 2 


DESCRIPTION OF EXPERIMENT 

The primary purpose of this experiment was to investigate the validity of the 

assumptions made by Eisenhauer and Spencer in deriving the curve of G (gj), defined 

s 

as a scattering function such that D°° = <*(X e ) ^G g (uJ u ) + G^u;^) , in which co u is that 

fraction of the upper hemisphere not subtended by the scattering surface, and is 
that fraction of the lower hemisphere not subtended by the scattering surface. This 
requires experimental knowledge of the variation of dose rate with solid angle with 
minimal structural or radiation complications and of the effect of wall thickness and 
structure shape on the dose rate. 


OUTLINE OF EXPERIMENT 

To achieve a direct dependence on solid angle, the detector was placed on the 
axis of a steel cylinder subtending solid angle u> = 1 - cos O, where G is the angle 
between the cylinder axis and a line from the detector to the upper (o? u ) or lower 
(o; u> ) edge of the cylinder (Figure 1). The dose rate could then be measured at 
different detector positions along the axis, expressed in terms of gj, and a com¬ 
parison made with G (gj). How'ever, the detector had to be placed below the level 

8 

of the contaminated plane to eliminate direct radiation, giving an expression 
D°° = a ^G g (uJ u ) - G g (c»J u ,) . The curve of G g (o;) was derived from Spencer's 
S a (d,oj) by the relation G g = 0.5 j^l - S a (d,aj)j. Since the G g (u>) curve in the Engi¬ 
neering Manual refers to fallout radiation, calculations were based on the S (d, g>) 
curve in Spencer’s monograph for Co-60 radiation using the formula 


D" - 0.5 Of [s a (u u ,) - S a (u u )j. 


Possible dependence of G (tu) on wall thickness was investigated by varying the mass 

s 

thickness of the cylinder; its dependence on structure shape was determined by re¬ 
placing the cylinder with a rectangular wall. 
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Figure 1. Solid Angles Subtended by Upper and Lower Limits of the 
Cylinder at Detector Position D 

However, the majority of the experiments were performed on the cylinder to 
take advantage of its symmetry, which permitted a point source to represent a ring 
source in the following way. If a source qS Q is placed a distance r from a detector, 
the response will be 



If qS is considered a linear radiation density, then the total response at the center 
° 2 

of a circle of such line charge is 2ETqS Q /r , since the contribution from one ele¬ 
ment is equal to that from any other. If this circumferential source is thought to 
represent an annular ring of area 2ardr, then the ring contributes 


dD = qS 

o 
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or 


dD 

dr 


2*tD 

m 


and the total field can be approximated by a summation of elements 2xrArD m , being 

represented by individual point sources a distance Ar apart along one radial line. 

Since 2-rrD = D is the basic datum, all measurements are expressed in this 
m exp 2 

form and are normalized to a field strength of 1 curie/ft . 

The same philosophy was followed in simulating the field for the rectangular wall 
except that the response would vary with the angle <f> between the perpendicular to 
the plate and the source-detector line. The total response would then involve a sum¬ 
mation over radius followed by a summation over the angle <f>. 


CHOICE OF MATERIAL 
3 

As Eisenhauer points out, his curve for G (w) is based on skyshine; i.e., radi- 

s 

ation scattered by an atmospheric medium. Practically, a denser scattering mate¬ 
rial was experimentally necessary and generally desirable for shelter applications. 
The significant choices were iron and concrete, both of which are presumed to be 
reasonably equivalent to air as scatterers, on a mass thickness basis, since the 
absorption coefficients of the lighter elements are quite close together at the photon 
energies of Co-60 or 1.12-hr fallout (around 1 MeV). Concrete would have greater 
applicability in structural shielding, but iron is easier to shape into the desired form 
and, more important, is much freer of voids and inhomogeneities. Since Monte 

4 

Carlo calculations indicate there is little actual difference in the scattering response 
of the two materials (Table 1), iron was chosen for the structural material. 


EXPERIMENTAL STRUCTURES 

The basic cylinder was constructed of steel l/2-in. thick, 5-ft high, and 2-ft in 
diameter. The height seemed to approximate an infinite cylinder closely enough, 
since even at the highest detector position was quite close to 0 (specifically 0.0188). 
This assumption was justified by varying the cylinder height and by observing the 
change of dose rate. The dose rate was close to its asymptotic value at a 5 -ft height 
(Figure 2). Thickness could be increased by adding 1/2-in. -thick half cylinders. 
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TABLE 1 


ROOF REDUCTION FACTORS COMPUTED BY THE MONTE CARLO 
METHOD FOR A 1. 25 MeV PLANE ISOTROPIC SOURCE ON A BARRIER 
OF 30 g/cm2 = 61. 5 psf EFFECTIVE MASS THICKNESS 


03 

Concrete 

Iron 

0.094 

I 

0. 0073 

0.0076 

0. 234 

0. 0170 

0.0176 

0.357 

0.0245 

0. 0257 

0.500 

0. 0316 

0. 0326 

0.658 

0.0360 

0. 0383 

0. 826 

0. 0391 

0. 0412 

0.913 

0. 0418 

0. 0424 

1.000 

0.0419 

0. 0429 



Figure 2. Variation of Detector 
Response with Cylinder Height 
for Various Source Positions 
(h = 3-in. depth) 
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shown in Figure 3, to the side of the cylin¬ 
der facing the source. Experiments were 
conducted for four cases of effective mass 
thickness: (1) 1/2 in. or 18.6 psf; (2) 

1 in. or 37. 2 psf; (3) 1-1/2 in. or 
55. 9 psf; (4) 2 in. or 74.5 psf. 

The rectangular wall was a 1-in. - 
thick steel slab 12 ft long by 6 ft high. At 
a detector position on a line 4 in. behind 
the wall, these dimensions gave solid 
angle subtended by the upper edge of the 
same order as that subtended in the cylin¬ 
der. The assumption that this would give 
a response close to that of an infinite plate 
was tested and verified by adding a small 
extra plate (Figures 4 and 5). 

The cylinder was placed over a hole 
of the same diameter and deep enough for 

■ * s s * c 
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Figure 3. 1/2-in. Thick Cylinders and Semi-Cylinders Used in the Experiments 


a 4-fit detector depth. The upper rim of the hole was lined with lead bricks to 
reduce edge effects. The rectangular wall was alongside a hole that only extended 
3 ft on either side of the detector position, thus blocking some of the wall from the 
detector line-of-sight. The change this made in the solid angle of the wall at the 
detector was significant only for the deepest three positions (Table 2). The rim 
was lined with lead bricks, and a steel plate covered the portion of the hole behind 
the detector to reduce airscatter. 
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Figure 4. 1-in. Thick Rectangular Wall with Dosimeters 

in Position 



Figure 5. 1-in. Thick Rectan¬ 

gular Wall Showing Detector 
Pit and Ratemeter 









TABLE 2 

CHANGE IN SOLID ANGLE DUE TO LIMITED DETECTOR HOLE 


Detector Depth 

L 

(in.) 

W 1 

W 2 

^1 

W 2 

36 

0. 0157 

0.0181 

0. 867 

24 

0. 0312 

0. 0335 

0.931 

12 

0.0787 

0. 0805 

0.978 

6 

0.1635 

0.1637 

0.999 


12 ft 



uj i = solid angle fraction subtended by the portion of the plate 
not hidden by the edges of the detector pit 

(i ?2 = solid angle fraction subtended by the entire plate. 
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POINT SOURCES 


Point sources were placed manually at selected radial positions. Sources of 
approximately 20. 7, and 0.5 curies were used, depending on detector depth and 
source radius. The 0.5 curie source was rod-mounted and positioned by hand; the 
7-curie and 20-curie sources were cranked into position through polyethylene tubing 
by a cable-type retrieval unit used for source positioning. 

INSTRUMENTATION 

At the deeper positions with smaller dose rates, the detector used was an ion¬ 
ization rate chamber removed from a standard Victoreen Model 592 Survey Meter. 

3 

This chamber, with an approximate volume of 320 cm , was modified by increasing 
the resistance across the chamber by a factor of 6 to provide greater sensitivity. 
The chamber was enclosed in a polyethylene cover, 1/8-in. thick, to ensure that 
electrons caused by free-air ionization did not penetrate the chamber wall. The 
output of the chamber was fed to a solid-state amplifier — a specially designed 
low-noise, feedback type. To remove statistical fluctuations, the output of this 
amplifier was fed through a time-integrating circuit to a John Fluke digital null volt¬ 
meter. In this manner, the output of the chamber could be read to approximately 
two parts in a thousand with excellent reproducibility. 

Near the surface, Victoreen Model 362 pocket dosimeters (200 mr) were used 
because of the appreciable vertical extent of the ratemeter. These were read with 
a portable charger-reader designed and constructed by Tech/Ops. 

CALIBRATION 

The ionization rate chamber used in the experiment was calibrated by exposing 
it to a standard Co-60 source previously calibrated by the National Bureau of Stand¬ 
ards. Both chamber and source were placed on an essentially massless calibration 
bench approximately 9-1/4 ft above, and parallel to, the ground. Voltage readings 

of the chamber-amplifier output were taken at various source-to-detector distances 

2 

ranging from 4 ft to 11 ft. These data were corrected for air attenuation, R atten¬ 
uation, and reflection of radiation from the ground to the detector. The resulting 
value of sensitivity obtained by extrapolating these calibration data to zero radius 
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for the detector-amplifier combination was 273 V output/r/hr. The detector re¬ 
sponse was determined to be essentially linear within 1/2 of 1% for voltage below 
10 V. Thus, when taking data for this series of experiments,we vary the source 
strength to keep the maximum detector output below this value as much as possible. 

A secondary calibration of the detector was performed by placing it within the 
hole and locating a small source of Co-60 in a fixed position relative to the detector 
immediately after the fundamental calibration was performed. This method of sec¬ 
ondary calibration was repeated to evaluate the detector performance at the start 
and conclusion of each day’s experimentation. 

The calibration of dosimeters necessary to convert microampere readings from 
the charger-reader to dose values in milliroentgens was conducted outdoors to re¬ 
duce scattering into the dosimeters. The dosimeters were carefully placed at heights 
above the ground that would minimize local effects of ground scattering without un¬ 
duly complicating the procedure. All exposures were to Co-60 sources of known 
strength. Groups of identical dosimeters were placed a given distance from the 
source. Different exposures were obtained by varying the time of exposure. The 
dose was calculated based cm an inverse square law behavior for source-dosimeter 
distance and no scattering, since the distances involved were small. 

OPERATIONAL PROCEDURES 

Personnel engaged in the experiment were equipped with 200 mr direct-reading 
dosimeters and film badges, and both the air-supported experimentation building 
and the control office were continuously monitored by Tech/Ops Model 492 
"Gammalarms. * In addition, during all manipulation of sources, two portable sur¬ 
vey meters were used for personnel monitoring. 

Dosimeters were accurately positioned at depths of 1, 2, 3, 4.5 , and 6 in. cm 
the axis of the cylinder. They were clipped to a plastic holder and exposed simul¬ 
taneously. While one operator monitored from a safe distance, the other attached 
the point source to a 14-ft handling rod and positioned it, in the case of the 0.5-curie 
source, or fastened the end of the polyethylene tube at the proper position and 
cranked out the larger sources. After the exposure was completed, the source was 
replaced in its storage container, and the dosimeters read and replaced for the next 
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run. Readings were taken at radial positions of 1.75, 2. 25, 2. 75, 3. 25, 3. 75, 

4.25 , 4.75 , 5.25 , 5.75 , 6.5 , 7.5 , 8.5, 11.0, 13.0, 15.0, 17.0, 19.0 , 22.5 , 27.5, 
32.5, 37.5, and 42.5 ft, and extrapolations were made at 1.0, 9.5, and 45 ft. 

The ratemeter was mounted on an aluminum rod and positioned by a pulley 
system. The same set of radial positions was measured at depths of 2. 75 , 4.25, 
5.75, 8.75, 11.75, 17.75, 23.75, 35.75, and 47. 75 in. The obscure depth settings 
resulted from an initial error of 0. 25 in. in the ratemeter position. The solid angle 
subtended by the detector varied with its depth as w = 1 - cos 0, where © is given 
by tan 0 = a/h; a is the radius of the cylinder; and h the detector depth. 

For the rectangular walls, dosimeters were placed at depths of 1.0, 2.0, 4. 0, 
and 6 in. and the ratemeters at depths of 4.0, 6. 0, 12.0, 24.0, and 36 in. The 
source was placed at distances of 1.75, 3.75, 7.5, 22.5, and 42.5 ft, each radial 
line being repeated at 0, 20, 40, 55, 70, and 80 deg. The angular positions were 
chosen on the basis of a crude estimation of flux change with angle. Solid angle 
measurements were based on the inner edges of the structures. 
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CHAPTER 3 


EXPERIMENTAL DATA 


Tables 3, 4, 5, and 6 show measured values of = 2 stD for effective 
mass thicknesses of 18.6, 37.2, 55.9, and 74.5 psf, respectively. Values at 
P = 1.00, 9.5, and 45 ft are extrapolated from curves of vs p. Table 7 gives 
D eX p for the rectangular wall. Dosimeter readings are at depths of 1 in., 2 in., 

3 in., and the first listed at (4a) 4 in. ; ratemeter readings are at 4(b), 6, 12, 24, 
and 36 in. Table 8 lists the value of solid angle fraction associated with each depth 
in the cylinder and wall. All values are normalized to a field strength of 1 curie/ft . 


TABLE 3 

DUSE RATE DATA DC 5-FT CTLDfDER 
(X, = IS. 6 psf) D.tp (tnr/hrVfcorWft) 


Xhfln.) 
P(ft) \ 

B 

B 

2.75 


4.25 

H 


6.0 

8.75 

1L 75 

17.75 

23.75 

35.75 

47.75 

1.00 

nr 

6.9 

5.2 

5.7 

4.0 

4.3 

2.79 

3.2 

L 87 

1.220 

0.696 

0.342 

0.149 

0.0790 

1.T5 

5. 2 

4. 4 

3.59 

3.65 

3.59 

2.94 

2.24 

2.29 

L 39 

0. 923 

0.464 

0.265 

0.111 

0.0667 

2.25 

3.9 

3.4 

2.92 

2.90 

2.92 

2.79 

L 70 

I L 71 

1.092 

0.741 

0.374 

0.217 

0.0940 

0. 0499 

2. 75 

3 2 

2.7 

2.33 

2.3 

2.33 

i L90 

L 44 

: L 40 

0.922 

0.621 

0. 320 

0.191 

0.0833 

0. 0439 

3.23 

2.9 

2.3 

2.02 

L 33 

2.02 

1.54 

1.26 

1.19 

0.907 

0.553 

0.290 

0 170 

0.0745 

0.0395 

3.75 

2.4 

2. 0 

L 73 

1.65 

1. 73 

1.39 

1.09 

1.06 

0.709 

0.490 

0.254 

0.149 

0.0690 

0.0365 

4.23 

2.1 

1.9 

L 55 

1.46 

1.55 

1.21 

0.972 

0.934 

0.643 

0.446 

0.231 

0.137 

0.0625 

0.0337 

4.75 

L • 

1.6 

L 32 

L 31 

L 32 

L 07 

0.944 

0.851 

0.563 

0.392 

0.204 

0. 125 

0.0674 

0.0317 

5.23 


L 4 

L 24 

L 22 

L 24 

0.996 

0.789 

0.769 

0.514 

0.360 

0. 199 

0.115 

0. C637 

0.0293 

5. 75 

L5 

L 3 

L 10 

1.07 

L 10 

0.873 

0.714 

0.689 

0.474 

0.326 

0.174 

0.107 

0.0506 

0.0274 

6.5 

L3 

LI 

0.997 

0.949 

0.997 

0. 732 

0.635 

0.611 

0.416 

0.279 

0. L.i 

0371 

0.0456 

0.0249 

7. 5 

L 2 

1.0 

0. 945 

0.954 

0. 945 

0.697 

0.554 

0.549 

0.371 

0.257 

0.139 

0.0967 

0.0409 

0.0222 

S.5 

LI 

0.9 

0.792 

0.756 

0. 792 

0.610 

0.497 

0.494 

0.332 


0.125 

0.0773 

0.0369 

0.0203 

9.5 

0.96 

0.9 

0.70 

0.66 

0. 70 

0.54 

0.439 

0.42 

0.2*5 

0.211 

0.117 

0.0702 

0.0339 

0.0190 

11.0 

0. 97 

0. 74 

0.579 

0.619 

0.579 

0. 494 

0.397 

0.412 

0. J6j 

0.196 

0.101 

0.0632 

0.030 

0.0166 

13.0 

0.65 

0.59 

0.515 

0.493 

0.515 

0.394 

0.327 

0.314 

0.217 

0.157 

0.0833 

0.0628 

0.025 

0.0139 

15.0 

0.60 

0.51 

0.435 

0.431 

0.435 

0. 349 

0.290 

0.278 

0.193 

o. m 

0.0744 

0.0469 


0.0124 

17.0 

0.53 

0.46 

0.399 

0. 595 

0. 398 

0.317 

0. 239 

0.253 

0.170 

0.121 

0.0661 

0.0410 

0.01*7 

o.oioa 

19.0 

0. 47 

0.40 

0.349 

0.336 

0. 349 

0.271 

0.225 

0.212 

0.152 

0.106 

0. 0598 

0.0365 

0.0177 

0.00988 

22.5 

0.39 

0.33 

0.301 

0.272 

0.3013 

0.223 

0. 192 

0.179 

0. 129 

0.0919 

0.0601 

0.0313 

0.0151 

0.00949 

27.5 

0.31 

0.27 

0.239 

0.226 

0. 239 

0. 178 

0. 149 ■ 

0.143 

0.103 

0.0729 

0.0407 

0.0255 

0.0121 

0.00704 

32.5 

0.25 

0. 22 

0.199 

0. 192 

0.199 

0.151 

0.129 

0.119 

0.0949 

0.0589 

0.0334 

0.0211 

0.0104 

0.00671 

37.5 

0. 23 

0. 194 

0.169 

0. 156' 

0. 169 

0.133 

0.1099 

0.1039 

0.0733 

0.0618 

0.0289 

0.0179 

0.0*694 

0.00604 

42.5 

0. 19 

0. 16 

0.142 

0.135 

0.142 

0.112 ; 

0.0914 

0.0899 

0.0649 

0.0467 

0.0263 

0.0163 

0.00784 

0.00442 

45.0 

0.19 

0.15 

0.133 

0.13 

0. 133 

0.106 

0.0927' 

0.064 

0.060 

0.0448 

0.0299 

0.0152 

0.00739 

0.00412 
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TABLE 7 


RATE DATA IN 6 FT z II IT RECTANGULAR WALL 
(X =3?.2p*f) D_(m/hrl/lcarie/ft) 


1.00 

0.28 

0.246 

0.228 

0.19S 

0.0720 

1.75 

0.211 

0.204 

0.164 

0.150 

0.0546 

3.75 

0.118 

0.119 

0.0813 

0.0677 

0.0241 

7.5 

0.0634 

0.0502 

0.0402 

0.0319 

0.0102 

22.5 

0.0135 ■ 

0.0144 j 

0.00876 

0.00860 

0.00337 

42.5 

0.0093 ] 

0.00887 ] 

0.00477 

0.00436 

0.00199 

45. 

0.0096 j 

0.0085 

0.00460 

0.00420 

0.00194 

h » 2 to. 


1.00 

1 

0.175 | 

0.169 | 

•.150 

0.151 

0.0650 

1.75 

0.134 | 

0.129 

O.lll 

0.112 

0.0451 

3.75 

0.0768 

0.0817 | 

0.0578 

0.0514 

0.0197 

7.5 

0.0410 

0.0412 

0.0289 

0.0235 

0.00809 

22.5 

0.0122 

0.0110 

0.00049 

0.00633 

0.00262 

42.5 

0.0070 

0.00629 | 

0.00347 

0.00336 

0.00173 

45. 

0.0068 

0.00610 

0.00339 

0.00332 

0.00166 


0.0375 

0.0305 

e.om 

0.00314 

0.00137 

0.00107 

0.00104 


0.0130 

0.0130 

0.0113 

0.00240 

0.00100 


1 

1.00 1 

0.1150 

0.105 

0.101 1 

0.104 

0.0548 

9.0120 

1.75 

0.0856 

0.0861 

0.0818 

0.0834 

0.0378 

9.0120 

3.75 

0.0528 

0.0537 

0.0403 

0.0387 

0.0169 

9.0113 

7.5 i 

0.0281 

0.0290 

0.0209 

0.0176 

0.00656 

9.90199 

22.5 

0.0089 

0.00759 | 

0.00409 

0.00531 

0.00202 

0.000985 

42.5 

0.0044 

0.00458 

0.00260 

0.0019S 

0.00120 

9.009734 

45. 

0.00425 

0.00450 

0.00252 

0.00182 

9.90116 

9.999722 


h = 4a in. 


1.99 

0.0850 

0.0760 

0.0750 

0.0800 

0.0280 

0.0130 

1.75 

0.0643 

0.0603 

0.0567 

0.0593 

0.0291 

0.0120 

3.75 

0.0360 

0.0401 

0.0301 

9.0296 

0.0137 

0.0109 

7.5 

0.0201 

0.0212 

0.0151 

0.0138 

0.00530 

0.90188 

22.5 

0.0072 

0.00591 

0.00395 

0.00405 

0.00169 

8.888695 

42.5 

0.0035 

0.00343 

0.00217 

0.#021« 

0.00196 

0.888598 

45. 

0.0034 

0.00332 

0.00208 

0.00206 

0.00104 

8.888558 





































TABLE 7 <ConS'd.) 

DOSE RATE DATA IN 6 FT x 12 FT RECTANGULAR WALL 
(X f - 37.2 paf) (nsr/hr(/<carif/ft) 


\ o 



| 

J,— 

T- 


6* \ 

0° 

i 20° 

1 

{ 40° 

j 

i 55° 

i 

70° 

J_ 

90° 

h = 4b m. 

s 

1.00 ! 

C.O&ZO 

| 0.09*0 

• 0.0*90 

; 0.0900 

| 0.0595 

0.0125 

1.75 

0.0734 

0.0752 

| 0.0651 

0.0633 

i 0.0351 

0.0125 

3.75 

0.03*1 

0.03%* 

0.0296 

0.0301 

0.0139 

0.0114 

7.5 

0.0202 

i 0.0206 

0.0251 

0.0139 

0.0054* 

0.00195 

22.5 

0.00650 

0.00577 

| 0.004+4 

0.00400 

j 0.00165 

0.00075 

42.5 

0.00365 

j 0.00363 

0.00230 

0.00220 

i 0.00096 

0.00057 

45. 

0.00340 

j 0.00350 

! 0.00221 

0.00212 

• 0.00092 

i _ 

0.00056 


b - 6 in. 


1.00 

0.0450 

T- 

‘ 0.0500 

0.0390 

| 0.0460 

i 

0.0250 

0.00900 

1.75 

0.0360 

! 0.0374 

0.0301 

0.0337 

! 0.0201 

0.00900 

3.75 

0.0212 

j 0.0215 

0.0177 

0.0196 

1 0.0101 

0.00999 

7.5 

0.0177 

0.0120 

0.00914 

0.00*53 

0.00356 

0.00122 

22.5 

0.00396 

0.003*2 

0.00296 

0.00363 

0.00110 

0.000504 

42.5 

0.00226 

0.00224 

0.00150 

0.00150 

0.00064 

0.00C396 

45. 

0.0021* 

0.0021% 

0.00145 

0.00145 

| 0.00061 

0.000359 


h - 12 tn. 


1.00 

0.0140 

0.0230 

J 0.0117 

0.0145 0.0120 

j 0.0105 

1.75 

0.01123 

I 0.0113 

0.009*5 

0.0116 0.00936 

! o.oioo 

3.75 

0.0076 

0.007*4 

0.0064* 

0.00721 0.00524 

0.00601 

7.5 

0.00476 

0.00497 

0.00376 

0.00370 0.00167 

] 0.000570 

22.5 

0.00179 

0.00142 

j C.000994 

0.000922 , 0.000462 

0.000182 

42.5 

0.00105 

0.00106 

0.0007S9 

' 0.0006*2 0.000292 

! 0.000163 

45. 

0.00103 

0.00103 

0.000750 

j 0.000679 ■ 0.000295 

J . L 

I 0.000160 


b = 203. 


1.00 

I 

| O.0O43O 

0.00440 

; 0.00340 

J 0.00355 

0.00359 

1 0.00339 

1.75 

0.00397 

0.0039S 

0.00314 

j 0.00340 

0.0037 

i 0.00320 

3.75 

j 0.00307 

: 0.00316 

0.0024.5 

; 0.00307 

0.00218 

; 0.00219 

7.5 

f 0.00201 

| 0.00204 

0.00156 

0.00157 

0.00069 

0.00025 

22.5 

0.000*2 

0.000*4 

0.00061 

0.00053 

0.00021 

' 0.00010 

42.5 

0.00050 ! 

0.00051 

! 0.00035 

0.000315 ! 

0.00013 

j 0.0000*5 

45. 

0.00049 

0.000499 

! 0.00034 

0.000309 

.. . i 

0.000125 

j O.OOOC79 

_i.. 


b = 3$ m. 


1.00 

0.00249 | 

0.00240 

1 "7- “T 

j 0.00139 i 0.00155 j 

0.00125 

| 0.00195 

1.75 

0.00211 

0.00214 

0.00135 j 0.00149 

0.00126 

1 0.00200 

3.75 

v.00145 j 

0.00150 

j 0.00126 j 0.00342 

0.0013S 

0.00209 

7.5 i 

0.00110 

0.00110 

1 0.000930 | 0.000912 ! 

0.000335 | 

0.000123 

22.5 

0.000465 

0.005495 

i 0.000336 ] 0.000292 j 

0.000104 

0.0000517 

42.5 

0.000272 

0.000283 

J 0.000193 : 0.000157 ! 

0.000^552 ; 

0.0000411 

**• _L 

0.000265 : 

0.000290 

0.000196 j 0.000156 

—---i ‘ 

0.0000540 i 
i 

0.0000402 
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TABLE 8 


SOLID ANGLE FRACTIONS, 01 , AND DETECTOR DEPTHS, h 
5 FT CYLINDER AND 6 FT x 12 FT STRUCTURE 



Cylinder 

Rectangular Structure 

h(in.) 

01 

^ / 

01 

oi , 


u 

u 

u 

u 

1.00 

0.0188 

ESI 

0.01426 

0.808 

2.00 

0.01821 

mm 

0.01388 

0.670 

2.75 

0.01777 

0.776 



3.00 

0.01767 

0.757 

0.01356 

0.556 

4.00 



0.01330 

0.465 

4.25 

0.01704 

0.666 



4.50 

0.01686 

0.649 



5.75 

0.01633 

0.568 



6.00 

0.01612 

0.553 

0.01250 

BESS 

8.75 

0.01497 

0.410 



11.75 

0.01367 

0.286 



12.00 



0.01120 

0.168 

17.75 

0.01165 

0.171 



23.75 

0.01007 

0.107 



24.00 



0.00862 

0.0710 

35.75 

0.00767 

0.052 



36.00 



0.00714 

0.0385 

47.75 

0.00610 

0.030 





















CHAPTER 4 


ANALYSIS OF EXPERIMENTAL DATA 
INTRODUCTION 

The calculational technique based on Spencer's data has proved very useful in 
evaluating and predicting protection factors of various shelters within the rather 
wide experimental limits and generally under conditions requiring uncertain exten¬ 
sion of the theory. It is therefore desirable to have an unambiguous test of the basic 
data and assumptions that have gone into the curves of the Engineering Manual. 

2 

It will first be necessary to outline the methods of the Manual, and their deri¬ 
vation from Spencer's work, ?nd then to demonstrate how the data that have been 
collected can be compared with this information. The method of calculating the 
infinite field dose rate from the experimental data and from the theoretical curves 
is described. The experimental formulas are given and applied to the data. 

DETERMINATION OF DOSE RATE AS A FUNCTION OF SOLID ANGLE 

The radiation reaching a detector located in a structure can be divided into 
three components: direct, air-scattered, and wall-scattered. The direct radiation 
is the straight-line component from source to detector; the component of air-scattered 
radiation that has been scattered through large angles is only slightly penetrating 
and of little consequence in a structure with thick walls and limited open ceiling; the 
small-angle component, introduced into calculations by the buildup factor, is near 
the primary in energy and penetration. 

The wall-scattered radiation has been deflected, by one or more Compton scat¬ 
terings , to the detector position and can usefully be considered as originating at the 
walls. Such an approach is valuable because the dose received at a detector can then 
be considered proportional to the solid angle of the walls viewed by the detector if 
the walls have a uniform distribution of radiation. Such an assumption is justified 
by assuming that the differences in distribution of impinging radiation, above and 
below the detector line, will be masked by multiple scattering in the wall; therefore, 
equal solid angles subtended above and below the detector line will yield equal detec¬ 
tor responses. The Engineering Manual curves are based on this assumption, which 
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would be expected to introduce the greatest discrepancy between experiment and 
theory. 


If this assumption is accepted, the dose rate at the detector in a given structure 

will be a function only of the solid angle subtended by the upper and lower limits of 

the structure. Then a scattering function G (oj) can be postulated such that the dose 

s 


rate D* = a(X e ) ^G s (a> u ) + G g (a> u/ )], 
and shape of the structure. 


where a is a function of the mass thickness 


3 

Eisenhauer has argued that wall-scattered radiation is close enough to air- 

scattered radiation in angular distribution to be approximated reasonably well by the 

already calculated distribution function for air-scattered radiation. Then the wall- 

scattering function G (w) is defined as 

8 


C S M = 


V 

I 


1(3', cos 6) d(cos 0) 


-(1-oj) 


(1) 


£(X, cos 0) is plotted on p. 84 of Spencer's monograph. 1 More directly, the curve 
for G (cd) found in the Engineering Manual is derived from that for S ( 3',w) in the 

S 8 

Monograph by the relation 


G s (cj) = 0.5^1- 

The structural constant a has been broken down into B w (X e )S w (X e ) E. A barrier 
shielding factor for the attenuation of radiation by the wall is B^. The curve in the 
Engineering Manual is the same as that for W(X,d) in the Monograph multiplied by 2. 
A weighting factor S w is described on p. 32 of Eisenhauer's explanatory treatise 
and represents the fraction of radiation scattered by the wall. Chart 7 of the Engi¬ 
neering Manual plots this weighting factor. 

The shape factor E is described by Eisenhauer on p. 31 and is plotted in Chart 
8 of the Manual. This factor varies from 1 for a structure with infinitely long walls 
to \ f~2 for a square structure. The variation can be interpreted as the additional 


«.( 3',u>) 




( 2 ) 
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field that would be viewed by a detector in a less eccentric structure through the two 
wall8 (in a rectangle) that would be removed to infinity in the structure with walls of 
infinite length. 

Thus a complete system of evaluation has been worked out based on the calcu¬ 
lated distribution of radiation intensities scattered into a particular solid angle by a 
scattering interface. This is expressed in the curve for G (u>) or S ( 3',w) modified 
by expressions for wall thickness and structure shape. 

COMPARISON OF THEORY WITH EXPERIMENT 

If these theoretical calculations have basis in fact, experimental knowledge of 
the variation of dose rate with solid angle should allow the comparison of those values 
with the infinite field dose rate computed from the curve of S (3’,u;). The simplest 
method of comparison is to insert the calculated values of o; u and into the Engi¬ 
neering Manual curve and to compare the result with the measured value. This gives 
an estimate of the reliability of the composite expression 

D°° = 0.5 S B E |S (cj , ) - S (u> ) 1. 

w w ^ a' u ' a' u'J 

Direct testing of the theoretical values of S (uj) requires additional assumptions 

3 

to derive an experimental curve. A discussion of this direct testing can be found in 
Chapter 6 of this report. 

CALCULATION OF INFINITE DOSE RATE 

The ideal situation involves infinite structures surrounded by infinite fields of 
radiation. The field is in actuality limited by time and the bounds of the testing area. 
Experimental determination of a rapid convergence toward complete response with 
structure height, expressible by the approximation of the solid angle subtended by 
the upper limit to zero, indicated that the practical limits of the structure were 
adequate. 

The problem of a finite cylinder can be accounted for in the comparison with 
S (3',w), but the contribution of the field beyond 45 ft amounted to a third or more 
of the total dose rate and required careful estimation based on theory and on data 
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from the field within 45 ft. As has been discussed in Chapter 2 of this report, a 
plane circular field about a cylinder can be represented by a succession of radial 
points. The measured dose rate from each point is multiplied by 2*rAr to give the 
dose rate from an annular ring. The total dose rate from a plane source is then 
found by summing these rings. The only inaccuracy lies in the finite size of Ar. 

An additional inaccuracy is introduced in the experiment with the rectangular wall 
when the field is simulated by a succession of radial lines at different angles to the 
wall and is separated by a finite A0. 

As was shown, the contribution from each annular ring is equivalent to the deriv¬ 
ative of the total dose with distance dD/dr at that radial distance. Then a point 
source qS Q at a distance r from a detector on the axis of a cylinder will represent 
an area source to give a gradient 


dD 

dr 


2*qS o e‘* lr B<Mr) r 


[ g s (w u> * g s (t V > 


( 3 ) 


where g g (a>) is a general scattering function, assumed independent of r for large r, 
and B(pr) is a buildup factor generally dependent on r. 

If a functional form can be found for B(pr), Eq. (3) can in theory be integrated 

over a range of r through which r is large enough to have given Ag its independence 

8 

of r. Then the contribution from this area to the dose rate will be 


D < r r r 2> 


-f 


2rqS B(pr) 

--- e^ r dr • Ag 


( 4 ) 


If it is found that Ag g reaches its asymptotic value at a distance equal to, or less 
than, the greatest experimental distance r , then r 1 can be set equal to r ,r_ = °° 
and 


D 


OO 


= D(r = r Q ) + 25rqS Q 


OO 

I 


e' Mr 

B(pr)dr • Ag g 


r 


o 


( 5 ) 
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where D(r = r ) is the dose rate from a field out to r . Thus the infinite field dose 
o o 

rate can be calculated if B(pr) is known analytically and is reasonably simple and 

if the asymptotic value of Ag can be found. 

s 

In traversing a medium, in this case air, a beam of radiation is attenuated by 
two principal mechanisms at this energy: (1) photoelectric absorption, in which the 
photon is permanently removed from the beam, and (2) Compton scattering, which 
can occur a number of times and which may eventually scatter the photon back into 
the solid angle viewed by the detector. In a narrow-beam attenuation experiment, 
the total absorption coefficient for a particular material and photon energy is meas¬ 
ured, but this coefficient in fact is made up of an energy absorption and a scattering 
coefficient. The energy absorption process follows a straight exponential decrease 
for each element of the radiation field, but this must be modified by a factor repre¬ 
senting the radiation scattered back into the beam, the buildup factor. Thus, 


D = D e 
o 




2 ' 
r 


5 

Experiments on the air-scatter of Co-60 radiation have led to a least-squares, 

0 72 

best-fit linear approximation of the buildup factor B(pr) = 1+0. 7(pr) * . This 

was chosen as one expression for B$*r) with p =p^.. As another extreme, B(pr) = 1 
was chosen, assuming that no scattered radiation rejoined the beam; a third choice 
was B(pr) = 1, p = p. , assuming no scattering took place at all. It was hoped that 

OO 

the experimental values of D found by each of these approaches would be close 
together and would bracket the correct expression among them. The three approaches 
then were: 

0 72 

1. B(pr) = 1+0.7(pr) ' p = p^. = total (narrow beam) absorption coefficient 

2. B(pr) = 1 P = P T 

3. B(pr) = 1 p = p A = energy absorption coefficient. 

As a derivation of asymptotic scattering function g (w), Eq. (5) can thus be integrated 

s 

and used if Ag g can be calculated from the data available and is found essentially in¬ 
dependent of r at the furthest measurement made. 
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The dose rate at a detector on the cylinder axis from wall-scattered radiation 
may be written: 


oj 


u 


2t 


dD 

dr 


2mjS 0 I d “ I 

to . 0 




( 6 ) 


where 


o 

q 


= radial distance from axis to source 

= RHF = 14.0 r/hr/curie at 1 ft from Co-60 

2 

= source density, curies/ft 

.2 


2*qS = 87.92 for q = 1 curie/ft 


o 

P = 
P = 


.-1 


absorption coefficient, ft 
slant distance from the source to the scattering area, ft 




2 2 a 2 a 1 1/2 

r + n +-g- 2 sin q (h cos 0 + r sin 0 cos <f> 

sm 0 8 


t> 


e 

h 


a = 


cylinder height 
detector depth 
cylinder radius 




GJ 


fit 


= 1 - 


= 1 - 


dose scattering function, steradian 
h + e 


-1 


J (h + e) 2 + 
h 

Jh 2 + a 2 
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As r becomes much larger than a and p ~ r, independent of G and 0, then 
Eq. (6) becomes 


dD 2»rqS o e-* r B(Mr) " f u a 2 * 


3F 


I cobQ § <x(X e ,0 0 ,Q,<t>) COS e du d<t>. 


( 7 ) 


V 0 


The integral in Eq. (7) is independent of r. If a function g (X , u>) is defined 

s c 


such that 


u 2r 


= 1 ds~e J «(X e .e o ,e,«co*e o dwd0, 


( 8 ) 


then 


<g _ 2>qS 0 e- |,r B(M r) p 
dr r 


*a< X e'"u>-8s< X e'V> 
L J 


(9) 


where g g (X e> c«;) represents the scattering properties of the part of the cylinder 
above cu, in a field of radiation beginning a distance r Q from the cylinder axis. 

The quantity actually measured in this experiment is dD/dr. If the assumptions 
made above are correct, the experimental data for a given detector position multi¬ 
plied by r/ should asymptotically approach a constant value for 

radii large with respect to the cylinder radius. This value is Ag , which can then 

8 

be used in Eq. (5) to find the infinite field dose. The result of this computation for 
the 1/2-in. thick cylinder is shown in Figure 6. The curves ire found to approach 
an asymptotic value as r increases. 

Table 9 lists Ag and D°° for two cylinder wall thicknesses, showing fairly close 
agreement among the three approaches and the expected bracketing of the dose rate 
incorporating buildup by the two extreme cases. 
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TABLE 9 


ASYMPTOTIC SCATTERING FUNCTION AggfX^u;) AND INFINITE-FIELD 
DOSE RATE FOR CYLINDER WALL THICKNESSES OF 18.6 AND 74.5 PSF 
USING DIFFERENT FUNCTIONS FOR B(fir) AT DETECTOR DEPTH OF 11. 75 IN. 


Scattering 

Function 

and 

Dose Rate 

x e (psD 

Bflir) = 1+0.7(per) 0,72 

M =^ T 

B#*r) = 1 

M =M t 

BQir) = 1 

Ag s 

ao 

18.6 

0. 0214 

0.0239 

0. 0230 

D 


11.17 

10. 26 

11.54 

Ag s 

74.5 

0. 00959 

0.0109 

0. 0105 

D°° 


5.13 

4. 75 

5.33 
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CALCULATION OF D°° BY THE ENGINEERING MANUAL METHOD 
The basic equation accounting for wall-scattered radiation is 


§- = B S EG (w )-G (u> ,) 
D w w s' u 7 s' u 7 


or 


k- = 0.5B S E S (o) ,) • S (u ) 

D w w a' u 7 a' u 7 

o L 

However, the experimental measurements include air-scattered radiation expressed 
as the buildup factor in the derivation of the far-field dose rate that must be taken 
into account as 

D r "I 0.088 B r "I 

irf ■ V 1 - S w> [><“„> - G a«V )J ■ <* - S w» [W > - S a<"u>} 

( 10 ) 

since there is no ceiling scatter contribution. The complete equation then becomes 


= 0.5B S E [s (u ,) - S (w )] 

D q w w ^ a' u 7 a' u 7 J 

+ 0.088B (1 - S ) | S (w ,) - S (w ) + 0.088 w . 

w' w 7 la' u 7 a' u 7 J u 


( 11 ) 


The last term represents the radiation scattered directly through the top opening of 
the cylinder. Since it is very small, less than 3% of the total even at the highest 
position in the thickest cylinder, it can be neglected. 

Equation (11) implies quite a difference in approach between the experimental 
and theoretical calculations for air-scatter effects. Experimentally, air-scattered 
radiation was handled as a multiplicative factor modifying the radiation flux at the 
scattering wall, ignoring air-scattered radiation striking the detector directly. The 
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Engineering Manual considers this radiation as an additive term dealing with direct 
air-scattered radiation only, attenuated but not scattered by the intervening wall. 
Since it would be difficult to separate experimentally the effects of air-scattered and 
direct wall-scattered radiation, the Engineering Manual approach doesn't exactly 
mirror the physical situation. Theoretical calculations indicate that skyshine is 
about 10% of the total for the thinnest (1/2 in.) walls; experiments indicate 6%. 

There is closer agreement for thicker structures, pointing to a maximum error of 
4% through accounting for skyshine. 

Equation (11) can be reduced to 

§“ = B w jyo> u , ) - S a (o> u )J [0.5S w E + 0. 088(1 - S w )J (12) 

for final calculations. 


PRESENTATION OF EXPERIMENTAL FORMULAS 

Dose rate was measured out to 42.5 ft, and an extrapolated value at 45 ft was 
used. As shown in Figure 6, Ag (X ,«) appears to have reached a constant value 

8 C 

90 

at this distance; also, the greater part of D (~ 2/3 of it) is contributed by the radi¬ 
ation within 45 ft. This is particularly important as the greatest source of error 

seems to be in the determination of Ag . 

s 

Another major problem, that of buildup factor, was handled by attempting to 
bracket it by three sets of assumptions: 


1 . 

V- 

- "total ~ X - 445 ft 

B(pr) = 1+0.70ir) 

2. 


^total 

BflJr) = l 

3. 


'"aba = 947.2# 

B(pr) = 1 


6 

where X is the mean free path in air. 
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This leads to three separate solutions of the integral in Eq. (3): 


*H*r] 

L 45 


-°. 2 8 -~ 

e dr = 220.7 


I ^ 


-r/445 


Idr = 2'qS 0 (^) = 159. 


n -r/947.2 x 

*■*. j —;— ldr - ^oCsso) = 221 - 3 ' 


and three expressions for D° 


d" = D(p =45) + 220. 7 |gl<u> u ) - g*< V ) j 
D” « Dtp = 45) - 159.4 [g 2 (u> u ) - g 2 ^,)j 
D" = D(p = 45) + 221.3 ^g 3 ( Uu ) - g 3 (lv ,j. 


Finally, Eq. (6) must be put into three forms: 


[**<«„> - g*<w u , >1 --——i- 

L J 87.92 e- r/44S |^1 ♦ 0. 7/445(Mr,°- 72 1 

[ g B<"u> * *»<“«' >] “ — '-r/445 3F 
[ g »'““ > ' g “' V ’] = 87.92 e^ 947 - 2 ^ 


to find the proper value of g (w,r — «). 
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Since the response of the wall would be symmetric about the perpendicular to the 
wall, measurements were made only in one quadrant and the resulting dose rate 
multiplied by four to simulate a structure with two infinitely long, parallel walls. 

The solid angle subtended by this structure was computed using the formula pre¬ 
sented on p. 68 of Spencer’s Monograph for one surface and multiplying this by two. 
The effect of the smaller size of the detector hole or the solid angle was found by 
treating the resultant complicated shape as a composite of simpler ones. 

EXPERIMENTAL AND CALCULATED VALUES OF D°° 

The raw data, D ex p(r), is presented in Table 4. In Table 10 D(r) is multiplied 
by At and is summed from 1 ft to 45 ft to give detector response from a plane source 
in that area;i.e., D(r = 45). In Table 11 D^^ • r/2rqS Q e ^ r B(pr) is presented for 
the 1/2 in.-thick cylinder, and is plotted against r in Figure 6. Tables 12 and 13 
present the asymptotic values of Ag , D(r = 45), and the values of D°°(oj) corre- 
spending to the three Eq. (15), for each value of X and structure shape. Calculated 

oo ® 

values of D /D q are presented in Table 14 and are compared with the experimental 
results in Table 15. 

The variation of D^cj) with u> in the 1 in. cylinder is shown in Figure 7. An 

OO 

extrapolated value of D at ui = 1 can be derived from this plot to give an experi- 

r 

mental estimate of B I 0. 

w [_ 

method appeared to be the approximation of the dose-rate points by an analytic func¬ 
tion of cjj as described on p. 45. The variation of dose rate with solid angle pre¬ 
dicted by the theory is compared with the experimental values in Figure 8. 


5 S E + 0.088(1 - S ) 
w ' w' 


however, a more consistent 
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TABLE 10 

D(r*45). SUMMATION OF DOSE RATES FROM ANNULAR RINGS 
OF RADIATION FROM 1 TO 45 FT; 1/2 IN. CYLINDER 





4 - 


CM 



f 


Q o 

II « 


32 


N 




U 


t C T T 


S 


= do*e rate at r; D2 = doee rate at r * 5r 





























TABLE 10 (Cant’d.) 

D(r=45). SUMMATION* OF DOSE RATES FROM ANNULAR RINGS 
OF RADIATION FROM 1 TO 45 FT; 1/2 IN. CYLINDER 


oooooooooooooooooooooooolc 
o’ © o © o o* o o o’ o o’ o’ o* © d o’ ©*©o© d o’ o d| d 

o> © ao ao co 

ficnnt'CHt'»ionNinnic«nHtfirtHi*nnt'N 
t-<Np**aoiQCM©ao©r)^*i-ioeo«“‘0©^-©©*r^‘ 

oooooooc5oooooooooooooooo 

dodo dodo o’ d d o d o' o’ dodo o* o* o d o 

iOMCO»OOOC)QOOOMONf<JDO«fO , 1 , ^fp0^aOp 

I- r -4 ««* C* <M © I- t£ l£> CO <D lO iQ t- CM I- (D CC> 00 »H © t- 

8SlS888c!oSSooSoSSoSS§SSo!S 

ddddddddoddddddddddddddd*'* 

I CO lO *■* 
lO©OWCO©lO^©CM©*f©ir5l~Ot^M'©<M©CO© 

©6J<*oo*-«inOiOMaftcocr>iOfMi>co^cc©,*9* H $fr5i'" 

COO<lOI-t-«rO»blP55CO<OC064C4M»H*-HH»HppO 


iOCMCOiO©OC>ao©©M©NraD©«f©«r*f©©r*aO© 
t- *-< ««* c* <m © I- t£ ir> CO ao to c~ iQ t- CM h- I- (D <£> ao »h © 

88lS888c!o8SooSoSSoS88SSo 

ddddddddoddddddddddddddd 

(0 io h 


•H fH O O 


I'lfififilPJtfOrtnNN 

OOOOOOOOOOOO 


o o o o o 


oooooooooooooocooooooooo 


o O © 00 © © 



§ © © 

h- © »H 

O fH ^4 


ao © © 

hl-®Nt0«y5« 

OPt'HtTrtftdOn 

©©.MiH^©©© 


0»H*^000*HfH^poo 

o' o o o o d o* d d o d o* 


ot-oaooaoo’i'COiorHtD 

^©^©COOOCOcDCOOt^lft 
OOOOOOOOOOOO 

d d o o' o o d o d o o d 

<N © © 00 h» 

•^A*«fOCOMar>onOioc>OJ<HiM«fi-.S'rfot^iopaoco 
c^ot-trjcocMoOpN’^'coc^irjcrin’^cvjooaooiftcop 

’«fCJ»'**"**H*-4fH©©FHFHF-4^li-4»-**HFM*HF^CMfH*-4F-t© 

ddddddddoddddddddddddddd 


lf5^rO<OMMCvj»Hi-4»HiM^^Hf-*000 


i 1 10 t ® N 

N N ift |> H 

8 8 1 8 8 


oooooooooooooooooooooo 


o 


' fH 

*1* 

© 

an 

CM © 

CD 

lO 

CM 

or 

CM 

■o* 

o < 

1 8 

8 

© 

CM 

00 

& 8 

8 

** 

CM 

CM 

M 

8 

■M* 

CO 

8 

8 ! 
o < 


o 


Nf © CO 


MfHt-MOO©0<OMoniOM©i^t-.oOtnpC'4if5iO©tO 

8 ooooM©«Hi-*ro©'<rcMC>t-^rcM'i©oD©©3 ‘ o* 

©10i0^'*C0WC0«CMCM«-*»-*fH»-4."* fH©©©©© 

o' d d o’ o' o o' o d o o o' o' o d o d d 


r*40t^'<t , coMcoonoo*rc^pp©*ot£>p 

nNN(onf)fO'rt’frt«’rii , ) , T«rt»J 


If) Ifl N V on If5 !' Cl 
©CMO0&I-©©© 
H H »4 O 


CONCOOOiOOJfHptCO ? S 8 
OiOr^t-’^OaDtO’^^Ot-Op 
COC0<OWMCM«-»»h*Hf"4©©©© 


O o o o o o 


to in io © »o t- © *h *r © cm © © cm irsiopt- 

S 0>OD'f©O*T©©QO<DCMir><NttC>CO©OO©iniO©** 

m* o* dodo o o' d o' o* o' o* o* o o* d d d o* o* d d o* 

-----—-—- - --off- 

M^*0<5f>©OFHt»tcCOplOCOtOO*H«-<*^C^ 

, tC>o©M©©i-«e'JvooNiftiM«$*y5co©«£cOfiOao 

C**© lift £j*H©CDO0t'-©U , > 1/5 *P , $ , C0&W<M*“*^4 fM*M»-«© 














































t—h/f 


a 

2 



*3 <a 

9 2 

ll ll n 

H H < 
a a. a 

ii u n 

a a a 


N 

o 

#*s 

U 

a. 

r- 

© 

+ 


ll 


li 


L* 

m 


a a 

fl PQ 


34 


1/947.2 ft. 





















TABLE 11 (Cont'd.) 



«j « N 

UJ l» 

. 3 a 

> ^ > 

it n 

H H 


§ 


a 

ii 

a 


n 


< 

a a 

ii n 

a a 


eg 


t- 

© 

+ 


Ii 

t: 

PQ 


II 


II 


w> u 

a. a 

PQ PQ 






35 







TABLE 11 {Coat'd.' 


f&V-h 0/39 



43 

i 

ii 

H 

a 

(I 


04 

fr- 

o 

t: 

*■» 

o 

♦ 


II 

H 

a 

ii 


43 

CJ 

t- 

s 


< 

a 

ii 


& 

o 


CQ 


>< 

a. 

a 


30 























1/2-in. CYLINDER 


CO 

O 

O 

a 

O 

O 

CO 

0 

0 

r- 

-f 

r-< 

O 

O 

0 

40 

O 

O 

*r 

0 

*4 

O 

O 

Cl 

t- 

0 

0 

<0 

1 - 

0 

0 

9* 

CO 

O 

O 

40 

00 

§ 

O 

O 

© 

$ 

O 

O 

O 

Cl 

s 

QO 

O 

Cl 

0 

0 

» n 

Cl 

0 

0 

CO 

*4 

Cl 

§ 

00 

*4 

CM 

O 

O 

SJ 

CM 

§ 

3 

r| 

O 

■> 

Cl 

0 

0 

«r 

CM 

CM 

§ 

S 3 

d 

s 

d 

O 

0 

d 

O 

O 

0 

0 

O 

d 

O 

O 

0 

0 

d 

0 

O 

O 

> 

0 

O 

0 


w cc t- t* i/i n t« Hh-nowwot'-ww f f ig 

HHHHHHHHHHHflNCIONMM I N N N 

OOOOOOOOOOOOOOOOOO •' o Q O 

ooooooooooooooooooc o © © 

OOOOOOOOOOOOOOOOOO c ooo 


P) (O f N tf 0^0 

ov to o p £ © p p 

H rt H N ' H fj fl 

OOOO J o o o 

OOOO J O O O 

dddodddd 


f-4^#-*©*-4*fcicJcoQ*H©©in©<*>pu?©««fp 
ci-rpr-pp^cico-riooaof'-oocoO© *?©©© 

oooooooooooooooooopoop 

OOOOOOOOOOOOOOOOOO OOOO 

OOOOOOOOOOOOOOOOOO OOOO 



— « 


o <Z> Cl Cl «-* O or; 15 

co m o co ■** © r f> © 

1.0 1.0 o o p p p co 

oppoooop 

00000000 


cigoco©o«i , «p-fcoco©©«H 

O vl m © *>h © Cl CM Cl ^ *f CO CO 

I-* h !• I« 00 <■/) « tf) « (13 (f) O) » 

oopoopooopopp 

0000000000000 


00000000000000000000 _ 


ssassasis 

i- i- r- i- 1 * r» h t- 

ooppppoo 

00000000 

© © o d © d o* o 


«s 

4S ** 

3 2 

fl I) 

H < 
a. a. 

H (I 


irt 10 1/5 10 m in in «o i.n 

t'* Cl I •* Cl f- Cl I- Cl t- 


t-0 m »0 © © © © © lO iO lO O tO 


h ci n n n t 


m in (O I* H o Ifj t- C N N N |» N 

HHHWHMNMCOf 


u u * 

5 a 4 . 

a » « 


0.7 (fir) 1 















T.U»L»: l-' 

ASYMPTOTIC SCATTERING FUNCTION AND INFINITE FIELD DOS’ HATE 
CYL1NDHICAL STIUCTCRES 
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J 53.48 
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28.69 
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! 15.c: 

j 0.088 

42.72 
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! 39.A9 
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! 0.072 

35.46 
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5 
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2A.0S 

0.045 

26.36 


4.5 

16.6A 

O.OA29 

26.11 

0.0A7 * 

24.12 

0.044 

26.38 


5.75 

12.A6 

0.0336 

19.88 

0.0375 ! 

18.44 

0.035 

20.21 | 


6.00 

13.02 

0.0352 

20.79 

0.037 ! 

18.92 

0.034 

20.54 


8.75 

7.623 

0.0210 

12.26 

0.023 ; 

11.29 

0.022 

12.49 J 


11.75 

4.660 

0.0138 

7.706 

0.0155 | 

7.131 

0.015 

7.979 


17.75 

2.267 

0.00701 

3.81A 

0.00795 | 

3.534 

0.0075 

3.927 


23.75 

1.302 

O.OOA20 

2.229 

O.OOA6 ) 

2.035 

0.0045 

2.298 


35.75 

0.5785 

0.00197 

1.013 

0.00225 i 

0.9372 

0.0022 j 

1.065 


47.75 

0.3102 

0.00105 

0.5419 

0.00120 j 0.5015 j 

0.00115 | 

0.5647 j 



38 


BUBLINCTON 


MASSACHUSCTTS 








tabu: u (i\.nt <3.) 

ASYMPTOTIC SCATTERING Kl‘NOTION AND INFINITE FIELD DOSE UaTE 
CYLINDRICAL STRCCTCRES 


i | "2 : -*• 

JL....J_ _ _i—_ l 


2 -in. Cylinder 


1.00 

21.70 

' 

... 

0.0537 

| 33.55 

0.05? 

| 30.95 

0.056 

34.09 

a. oo 

18.42 


0.0456 

1 28.48 

0.050 

\ 26.33 j 

0.049 

29.04 

2, 73 

16.04 


0.0421 

25.33 

0.047 

j 23.53 j 

0.044 

25.79 

3.00 

15.60 


0.0384 

! 24.07 

0.045 

22.77 ; 

0.042 

24.99 

4.25 

11.89 


0.03i 7 

« 1S.S9 

0.036 

17.63 

0.034 

19.41 

4.5 

11.78 


0.0295 

I 29 

0.035 

17.. o | 

0.032 

19.86 

5.75 

9.106 


0.0246 

14.54 

0.023 

* 13.57 ! 

0.029 

15.30 

6.00 

9.454 


0.0253 

* 15.04 

0.028 

5 13.92 j 

0.0265 

15.32 

9.75 

5.214 


0.0149 

j 8.502 

0.017 

j 7.924 j 

0.016 

8.755 

11.75 

3 2°** 


0,00959 

5.339 

0.0109 

i 4.959 j 

0.0103 

5.546 

17.75 

1.506 


0.00475 

1 2.554 

0.0054 

| 2.367 j 

0.0052 

2.657 

23.73 

0. So 23 


0.00274 

1.457 

0.0031 

i 13 < 7 1 

0.00295 

1.505 

35.75 

0.3722 


0.00126 

‘ 0.6503 

0.0014 

J 0.5934 | 

0.00135 

0.6710 

^ , 

0.2141 

| 0.0007*7 

! 0.387S 

0.00091 

1 - °^ 92 J 

0.00097 

0.4066 


TABLE 13 

INFINITE FIELD DOSE RATES 
RECTANGULAR STRUCTURE 


Depth 

(in.) 

oo 

D i 

1.0 

36.98 

2.0 

26.12 

3.0 

18.52 

4.0 

14.65 

6.0 

9.27 

12.0 

3.803 

24.0 

1.716 

36.0 

0.918 
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TABLE 14 

CALCULATION'S USING THE ENGINEERING MANUAL METHOD 
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TABLE 15 

COMPARISON OF EXPERIMENTAL AND THEORETICAL VALUES OF D 
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Figure 7. D°°(aj) vs a>, 1-in. Thick Cylinder 
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Figure 8. D* vs w, Calculated and 
Experimental Values 
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CHAPTER 5 


RESULTS AND CONCLUSIONS 
EXPERIMENTAL ERROR 

The opportunities for error lay in the interpretation and extrapolation of the data 
as well as in the experimental limitations of accuracy. In taking readings, error 
could arise from 

1. Incorrect source position 

2. Surface irregularities 

3. Incorrect detector position 

4. Consistent inaccuracy of the detectors 

5. Normal dosimeter spread and reader-charger readout limitations 

6. Finite size of detectors. 

Calculational errors would be most likely in determining the value of Ag from the 

s 

plots (see Figure 6) and in extrapolating D (w) to w = 1. 

Most of the experimental errors can be considered more or less statistical fluc¬ 
tuations because of the large number of independent positionings and readings tliat 
went into the final values of D°°. For the dosimeters this is true of 1, 3, and 5 above; 
for the ratemeter, which was left undisturbed while a radial series of points was 
measured, 3 could have been serious if accurate positioning of the ratemeter had not 
been ensured. Surface irregularities were recognized and compensated for. Since 
the same dosimeters were consistently used in the same positions, sizable differences 
in individual responses would distort the final result. However, these differences 
did not exceed 2 to 3% in subsequent measurements. The ratemeter, for which 6 
might be significant, was used only at the greater depths where the flux changed rea¬ 
sonably slowly. 

Scatter of individual measurements showed up strongly in the plot of Ag vs 

s 

radius in Figure 6 to give a somewhat uncertain asymptotic value. Hopefully, the 

use of three approaches, calculated from the Ag curves on different days in different 

s 
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moods, leads to some statistical reduction of error in the average value. The vari- 

OO 

ation of the quadratic fitted to the D (a;) points from those points was no more than 
2% at the larger values of u>. There could be up to 4% error in the treatment of 
air-scattered radiation. On the basis of these considerations, the data should be 
accurate to within 10%. 

COMPARISON OF EXPERIMENTAL AND THEORETICAL RESULTS 

As shown in Table 15 there is good agreement between theory and experiment 
for all but the thinnest wall and largest solid angle. In thin walls, of course, the 
assumption breaks down that multiple scattering will mask the difference in response 
to various parts of the wall. 

The G (cj) curve rests on the assumption that equal increments of solid angle 
s 

will give equal response no matter what the height of the angles above the source 
plane or equivalently the detector position below the plane. But in reality the 
expression 


r 

D = 0. 5 a I S (o> ,) - S (fa> ) 
la' u ' a' u 7 

is valid only for a detector position at the level of the source plane in a semi-infinite 
cylinder; i. e., 


D°°(u; u , =1) = 0. 5 or 1^1 — 0 = 0.5a. 

For positions below the surface, the correct expression would be 

D% u , ) = 0. 5 f(d)a [yu> u , ) - S a (a> u ) J , (16) 

with f(d) > 1, since lowering the detector is equivalent to raising the source plane 
thus putting it closer to each incremental scattering area that is still exposed. This 
form, of course, is necessitated by the retention of the theoretical S (u>) rather than 
by adopting a new, experimental form. The theoretical value would be expected to 
have sufficient generality to justify modification rather than replacement. In Chapter 
6 the experimental curves are found to tend toward Spencer’s as a limiting value. 
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This leads to two conclusions about the terms of Eq. (11): 

1. The values of the components of 0.5S B E+0. 088B (1 - S ) 

I w w w w 

may be wrong since the theoretical expression should be quite 
accurate for a> u , ~ 1. Then if the composite value is changed 
to agree with an extrapolated experimental value of D°°atoj = 1, 
a new set of calculated values for D(uj) is obtained as shown 
in Table 16. A comparison of <* eX p with the theoretical expres¬ 
sion is shown in Table 17. 


2. For each measured value of a?, f(d) can then be determined, 

assuming S (o>,d) actually represents wall-scattered radiation 
31 

and the wall-scattered and air-scattered (0.5 and 0.088) com¬ 
ponents vary identically with depth. Table 16 also shows f(d). 


CONCLUSIONS 


Agreement within about 20% is found for most of the positions and wall thick¬ 
nesses, except the thinnest, for which the greatest discrepancy would be expected. 


If the validity of S (d,uj) is accepted and the barrier factors are questioned, a 




different value for | 0. 5 S B r E + 0. 088 B (1 - S ) | is found from the known be- 

w w w w 


) i* 


havior of S at oj = 1 and the extrapolated value of D°°. The factor S (o>) is defined 
a a 

as equal to 1 at w = 1 and as equal to 0 at w = 0. Therefore, theory requires that a 
detector at the intersection of a semi-infinite cylinder with a plane of radiation, sub¬ 
tending lower solid angle w u , = 1 and upper angle w u = 0, receives dose rate 


D 


= or jyi) - S a (0) = a. 


Thus an experimental value for B j 0. 5 

L 

for the finite size of the cylinder and detector. The missing portion of the cylinder 
between 5 ft and infinity represents a very small part of the dose rate at the higher 
detector positions and can be adequately approximated by its theoretical value, 

QtS a( a5 u ) * 


S E 
w 


0. 088(1 - S w ) 


could be found, except 
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TABLE 16 

D® CALCULATED FROM EXTRAPOLATED VALUES OF D® TO DETERMINE f(d)* 

exp ' ' 


“V 

D°° 

calc 

D 

o 

D® 

exp 

D 

o 

f(d) 

D® 

calc 

D 

o 

D 

_ex£ 

D 

U 

f(d) 


l/2-in. Cylinder* 

2 

1-in. Cylinder 

0.917 

0.0947 

0.1076 

1. 136 

0.0943 

0.1092 

1. 158 

0. 836 

0. 0751 

0. 0917 

1. 221 

0. 0748 

0. 0895 

1. 196 

0.776 

0. 0650 

0.0795 

1. 223 

0.0647 

0. 0834 

1.289 

0. 757 

0.0598 

0.0766 

1.2S1 

0.0595 

0.0756 

1.270 

0. 666 

0.0482 

0.0628 

1.303 

0. 04SO 

0. 0668 

1.392 

0. 649 

0.0467 

0.0620 


0. 0465 

0.0594 

1.278 

0.568 

0. 0373 

0. 0501 

1.343 

0. 0371 

0.0515 

1.388 

0. 553 

0.0357 

0.0490 

1.373 

0.0356 

0.0483 

1.358 

0.410 

0. 0228 

0.0337 

1.478 

0.0227 

0. 0323 

1.424 

0. 286 

0.0149 

0.0236 

1.584 

0. 0148 

0.0215 

1.450 

0. 171 

0.00752 

0.01286 

1.710 

0.00749 

0.0114 

1.522 

0.107 

0.00407 

0.00787 

1.934 

0.00406 

0. 00651 

1.604 

0. 052 

0. 00175 

0.00372 

2. 126 

0.00174 

0. 00296 

1.697 


0.000856 

0.00206 

2.407 

0.000853 

0.00160 

1. 876 


’*** D%, - B w j^O. 5S w E • 0. OMO - S.> , - S^, ■ - a ) - S^J 

__ D *i0 , --P-1U , P .1U 

I - S <- ) 1-0. OteS 0. 9^2 

I IUD 


1. d“(I) - «2. M 

o = 62. M 


2 . D“,i, -«Lrr 

o « 62. 23 

- o. is: 
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TABLE 16 (Cont’d.) 

D* CALCULATED FROM EXTRAPOLATED VALUES OF D TO DETERMINE f(d) 

exp ' ' 


w u' 

oo 

calc 

D 

o 

D°° 

exp 

D 

o 

f(d) 

O 

o°g. 8 
o 

e>: 

exp 

D 

o 

f(d) 


1-1/2-in. Cylinder** 


2-in. Cylinder 

4 

0.917 

0.0812 

0.0930 

1. 146 

0. 0580 


0.0675 

1. 164 

0. 836 

0. 0644 

0.0774 

1. 203 

0.0459 


0. 0573 

1. 248 

0.776 

0. 0558 

0.0689 

1. 236 

0.0398 


0.0510 

1. 283 

0. 757 

0. 0513 

0.0674 

1.315 

0.0366 


0.0484 

1.321 

0.666 

0.0413 

0.0521 

1.262 

0.0295 


0. 0380 

1.286 

0.649 

0.0401 

0.0525 

1.311 

0. 0286 


0. 0368 

1.286 

0.568 

0. 0320 

0. 0400 

1.251 

0. 0229 


0. 0293 

1. 282 

0.553 

0. 0306 

0.0418 

1.368 

0.0218 


0.0303 

1.387 

0.410 

0.0195 

0.0247 

1. 264 

0. 0139 


0. 01711 

1. 228 

0. 286 

0.0127 

0.0155 

1.216 

0. 00910 


0. 01074 

1.179 

0.171 

0.00645 

0.00767 

1.190 

0. 00460 


0.00514 

1.117 

0. 107 

0. 00349 

0.00448 

1. 285 

0. 00249 


0.00293 

1. 176 

0.052 

0.00149 

0. 00204 

1.366 

0. 00107 


0. 001308 

1. 223 

0. 030 

0.000734 

0. 00109 

1.486 

0. 000524 


0.00078 

1.488 


1-in. Rectangular 5 

3. D*(l)-J 

ks.n 







a * S3. SS 


0. 808 

0. 0607 

0.0744 

1. 226 



15 


0.670 

0. 0420 

0.0526 

1. 252 

4. D*1) = M.*1 

a • 33.23 


0.556 

0. 0307 

0.0373 

1.215 



-j, » 0. CT*»2 


0. 465 

0.0229 

0.0295 

1. 288 

5. D*(l) - S5. *S 


0.340 

0.0151 

0.0187 

1.238 



157 


0.168 

0.00602 

0.00765 

1. 271 





0. 071 

0.00182 

0.00345 

1. 896 





0.0385 

0.000984 

0.00185 

1. 880 
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TABLE 17 


0.5S w B w E + 0. 088B w (l - S w ) : EXPERIMENTAL AND THEORETICAL VALUES 


X e (psf) 

a 

exp 

^theor. 

a exp 

Analytic 

Visual 

, a theor. 

18.6 

0.1255 

1 

0. 201 

0.624 

37.2 

0.1252 


0.177 

0. 707 

37.2 (rect.) 

0.1112 


0.119 

0.934 

55.9 

0.1078 

0.113 

0.129 

0.836 

74.5 

0. 0769 

0. 080 

0. 091 

0.845 


Because of the danger of direct beam contribution, the highest detector was an 
inch below the surface, subtending a solid angle = 0.917. This is sufficiently 
close to u , = 1 to permit extrapolation of the dose rate to that value by visual ex- 

CO QO 

tension of the D (cj) vs gj curve or by analytic curve fitting representing D by a 
quadratic in u determined by the method of least squares. These curves are com¬ 
pared with the data in Table 18. The results of the two methods are compared with 
the theoretical value in Table 17. A new set of calculated values of D°° is obtained 
that at least converges toward the experimental value at o> = 1, where agreement 
would theoretically be inevitable, though the discrepancy at greater depths is now 
larger. Values of f(d) for each thickness are found to form a fairly smooth family 
of curves (Figure 9). A possibly significant factor that has been neglected is the 
variation of S (d.tu) and S(d) with depth or detector distance from the standard height 

A 

of 3 ft. An attempt to estimate this has been included in Chapter 6. 

General agreement has, therefore, been found between the results of this experi¬ 
ment and calculations made using the method and curves of the Engineering Manual. 
However, significant discrepancies exist that can be accounted for theoretically by 
assuming that one or all of the factors B^, S w> and E are incorrect and that there 
is an additional factor, f(d) that varies with detector position and wall thickness. 
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CHAPTER 6 




S a (w): DERIVATION OF AN EXPERIMENTAL CURVE 

The infinite field dose rate has been found for positions below ground in a steel 
cylinder and behind a rectangular wall. These values have been compared with 
those predicted by the method of the Engineering Manual using the expression 

S- = o [s^.) - s a <u u )j. 

The factors B 0.5 S E + 0.088(1 - S ) are represented by a and can be deter- 
IV I W IV I 

oo 

mined either by using curves given in the Manual or by extrapolating D to its value 
at a> = 1 (see p. 45). 

Since the object of this experiment is to evaluate the calculated curve G (gj) or 
the equivalent S (w), it is desirable to eliminate the additional variables B , S , 
and E by following the second course mentioned above. This experimental value of 
a can then be used to find an experimental set of values of S (oj) through the equation 

A 


D . . 

D"< w u' ’ W u> 

S a (w u' > = - 2 “5- + W < l7 > 

if S (ci> ) can somehow be evaluated, 
a u 

For all values of « u , encountered in the experiment, u> ;j is small and so is the 
correction term S (w ) (corresponding to the difference between the experimental 

8 U 

structure and an ideal, infinite one). Further, this range of S a (o>) can be closely 
approximated by a straight line; i. e., S (o> u ) = bo^. Thus a simpler formula can 
be found with negligible introduction of error 

D®° 

D" ( "u' '“u> 

s,«V > - - 2 —;- + >>v < 18 > 


so 


I 


s 







The problem is now reduced to evaluating b. If we return to the original form 






» 


"u> = “ 


= a 


[ S a«V > 

[ W > 


W] 

b "u] 


and differentiate, we find 


dD 


oo 



- erb 


(19) 


( 20 ) 


since S (cj ,) is independent of cj . From the variation in dose rate,dD°°/doj can 
be found with a change in the height of the cylinder, corresponding to a change in 
o> u . At a given detector position, the height of the cylinder was increased in steps 
from 1 tc 7 ft, and measurements were made with the source at three radial posi¬ 
tions, using the 1/2-in. cylinder. By dividing the change in dose by the accompany¬ 
ing change in solid angle, values for dD(r)/do? u could be found for each radial position. 
An approximation of dD°°/dcJ u could be found by adding the three dD(r)/dw u • Ar and 
a fourth term representing the field from 45 ft to infinity (Table 19). This is obvi¬ 
ously a crude approach based on marginal data, but it is justified by the noncriti¬ 
cality of the exact value of b and by the consistency of the value calculated with the 

derived curve of S («). 

a 

The theoretical curve of S (w) used is taken from p. 120 of Spencer’s monograph 
and is based on a constant detector height of 3 ft. Since the experimental variation 
of u was achieved by changing the detector height in a structure of constant dimen¬ 
sions , a correction must be made to relate the experimental data to that which would 
be found at the standard height. 

Using Spencer’s definition, we find 


- 1 + u 

s a (“) = Sp) jj d(cos 0) /(d, 
-1 


cos 0). 


• U 


I 


c 
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The total air-scatter received from the upper hemisphere by an isotropic detector 
at height d is represented by S(d) and will decrease with height symmetrically about 
the source plane. The integral term gives the fraction of S(d) scattered into the 
solid angle of the upper limit and, of course, equals S(d) for gj = 1. Since detector 
positions below ground receive higher dose rates than positions subtending equal 
solid angle above the source plane, as on p. 45, it can be assumed that Spencer’s 
curve of £(d, cos 0) will continue to increase for values of d less than 3 ft. Because 
the value of £ seems to increase with height in the same proportion for all values of 
cos © less than 0, the assumption was made that the integral over £ will increase by 
the same amount. By plotting £(d - 1) vs d as given in Spencer’s curve (Figure 10), 
extrapolated values of £ for d less than 3 ft could be found and the corresponding 
values of 

-1+ w 

jj d(cos 0) £(d, cos 0) 

-1 

as shown in Table 20. Again, the inaccuracies of the approach are offset by the 
smallness of the correction and by the fact that any correction will be an improvement. 



Figure 10. I(d, cos 0 = -1) vs d 



TABLE 20 

S(d) and £( d, cos 0): CORRECTION FACTORS TO FIND D°° (3’, cos 0) 

s 


S ~ (3 '' u ’ )= f(f) X(i) S a< d, “) 


a 


d 

(in.) 

S(d) 

A(d) 

S(d) 

M3’) 

Correction 

S( 3') 

X(d) 

36 

mm 

m 

1. 00 


— 

-1 

1 

IB’- 

1.043 


1. 004 

-2 


0.429 

1. 043 


1. 002 

-2.75 

0. 088 

0.4295 

1.043 

0.959 

1. 001 

-3.0 

0.088 

0.4295 

1. 043 

0. 959 

1. 001 

-4.25 

0.088 

0.43 

1. 043 

0.958 

0. 999 

-4.5 

0. 088 

0.43 

1. 043 

0.958 

0.999 

-5.75 

0. 088 

0. 432 

1.043 

0. 954 

0. 995 

-6.0 

0. 088 

0.432 

1. 043 

0. 954 

0. 995 

-8.75 

0. 088 

0.4325 

1.043 

0.953 

0.993 

-11. 75 

0.088 

0.434 

1.043 

0. 949 

0.990 

-17. 75 

0. 0863 

0.437 

1. 023 

0. 943 

0. 963 

-23.75 

0. 0850 

0.44 

1. 007 

0. 936 

0.942 

-35.75 

0.0844 

0.448 


0. 920 

0. 920 

-47.75 

0.0835 

0.454 

0.989 

0. 907 

0. 898 


COMPARISON WITH THEORY 

Once the experimental values D°°(u;, d) have been changed to D°°(a),3), values of 

S (oj) can be found for each wall thickness (Table 21). The values derived from 
£L 

these data are plotted and are compared with the calculated curve in Figure 11. 
Curves through these data points are presented in Figure 12. Experiment and theory 
come into closer agreement as wall thickness increases, indicating that the major 
cause of the discrepancy lies in the assumption of multiple scattering in the inter¬ 
vening wall to give essentially isotropic scattered radiation. This assumption was 
made 3 to permit the use of the solid angle subtended by the wall at the detector as 
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TABLE 21 

S (u): EXPERIMENTAL AND CALCULATED VALUES 
a' 


C 0 

ess 

lin. 

1-in. Rect. 

1-1/2 in. 

2 in. 

Calculated 

0.917 

i 

IQ 


BQ 

m 

BQ 

0. 836 


ass 


bq 

H 


0. 808 



0.680 



OQ 

0.776 

0.643 

0.677 


BQ 

0.674 

hq 

0. 757 

0.620 

0.615 


91S 

0.639 


0.670 



0.482 




0. ‘'56 

0.502 

0.543 


0.492 

0.504 

■SB 

0.649 

0.502 

0.483 


0. 496 

0.488 

mm 

0.568 

0.406 

0.418 


0.379 

0.390 

■SB 

0.556 



0.343 




0.553 

0.398 

0. 394 


0.395 

0. 400 

BL 

0.465 



0. 273 



0.225 

0.410 

0. 275 

0. 265 


0.236 

0. 230 

0.185 

0. 340 



0.174 



0.145 

0.286 

0.194 

0.178 


0.150 

0.146 

0.122 

0.171 

0.105 

0. 0945 


0.0753 

0.0712 

0.0628 

0.168 



0.0745 



0.0620 

0.107 

0. 0648 

0.0549 



0. 0418 

0.0350 

0.0710 



0. 0343 



0.0218 

0.052 

0. 0317 

0. 0262 


0.0219 

0. 0201 

0.0159 

0. 0385 



0. 0195 



0.0112 

0.030 

0.0183 

0.0151 


0. 0127 

0. 0127 

0.0084 


the sole criterion of the radiation scattered to the detector, whether from above or 
below the detector plane. Obviously, this assumption will be more reasonable in a 
structure with thicker walls in which multiple scattering will more effectively con¬ 
ceal the different angular distributions of the radiation falling on different sections 
of the wall. 

























Figure 9, in which f(d) is plotted for each value of X e> gives an indication of 
the variation of the scattering properties of the wall with detector positions. This 
is essentially the same as a ratio of the experimental and calculated values of S (o>) 
given in Table 21. A series of curves for S (&>) has thus been found that, in the 

2L 

limit of great thickness, approaches that calculated by Spencer. They differ sub¬ 
stantially for the smallest solid angles and thinnest walls, though this discrepancy 
is somewhat offset in dose calculations by a discrepancy in a in the opposite direc¬ 
tion. Since experiment predicts a higher dose rate than theory, substantiating ex¬ 
periments would be desirable. 

This experiment indicates that the calculational methods of the Engineering 
Manual are entirely adequate for the evaluation of most shelter situations. However, 
there does seem to be room for a considerable refinement of theory to interpret the 
discrepancies encountered. 
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